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'^ Abstract 

a\ 

CN We review the role of zero-temperature entropy in several closely- 

^ related contexts in QCD. The first is entropy associated with disor- 

r^ dered condensates, including (G^^). The second is effective vacuum 

Y^ entropy arising from QCD solitons such as center vortices, yielding 

Oh confinement and chiral symmetry breaking. The third is entangle- 

(— I ment entropy, which is entropy associated with a pure state, such as 

'— ' the QCD vacuum, when the state is partially unobserved and un- 

,__( known. Typically, entanglement entropy of an unobserved 3-volume 

I> scales not with the volume but with the area of its bounding surface. 

^^ The fourth manifestation of entropy in QCD is the configurational 

\^ entropy of light-particle world lines and flux tubes; we argue that this 

^O entropy is critical for understanding how confinement produces chiral 

(^ symmetry breakdown, as manifested by a dynamically-massive quark, 

"^ a massless pion, and a (qq) condensate. 

• ^ 1 Introduction 

X 

j^ Entropy is a precise tool for specifying imprecision. In statistical mechanics 

and in certain areas of field theory, the entropy S{qi), a function of a small 
number of state variables qi, is defined as the logarithm of the number A/" ^ 1 
of system states consistent with the small number of state variables: S = 
In A/". The imprecision arises because all the system states have essentially 
the same values of the state variables, such as (in thermodynamics) total 
internal energy U and particle number A^, and one does not know or care 
precisely which system state one is dealing with. In field theories with soliton 
condensates, such as QCD, it also happens that there is a vast number of 



disordered configurations of about tlie same action contributing to tlie unique 
vacuum state, and it becomes useful, witli a sligfit abuse of notation, to speak 
of the entropy of the vacuum as a measure of the number of configurations 
of these disordered states. 

QCD, hke any many-body system, has entropy at finite temperature, 
where there are a great many pure states that contribute to the mixed ther- 
mal state. Although there are non-trivial problems in evaluating QCD finite- 
temperature entropy, having to do with infrared singularities and strong cou- 
pling, we will not discuss it further except to mention an analogy between 
thermodynamic entropy and the kinds of zero-temperature QCD entropy of 
interest to us. The canonical partition function Zj^ for an isolated thermal 
system at fixed particle number A^ is the trace over all A^-particle states: 

Z^ = Tre-^^ (1) 

Here /3 = 1/T is the inverse temperature and the Hamiltonian is H. This 
partition function can also be written as an integral over energy, with an 
entropic weight factor, yielding the Helmholtz free energy: 

Zr, = fdE e^(^'^)-/'^ ^ e-'^^^-^^). (2) 

The thermodynamic quantity U is the energy E at the sharp maximum of 
the integrand, where dS/dU = 1/T. 

Throughout this paper we mostly use the dilute gas approximation (DGA) 
to describe soliton condensates of QCD. In the A^-soliton sector of the DGA 
grand partition function Z the action of A^ solitons is the sum of A^ single- 
soliton actions Ic plus a quadratic term whose path integral yields a product 
V of determinants times integrals over the collective coordinates qi arising 
from zero modes — derivatives of the 0cj with respective to the qt — that 
produce zero eigenvalues in some of the determinants. Then 



cc 



where the sum over A^ is the sum over distinct solitons, cc indicates the sum 
over zero modes, and we suppress writing the determinantal factors. Note 
that all solitons have the same action, which needs no index i. The terms 
(1/A^!) '^^^'D can be regrouped into an entropy factor, much as in Eq. ([2|. 
In gauge theories such as QCD the action is ~ l/g"^ so that condensates form 



when the couphng g^ is large enough. There should be no non-perturbative 
condensates for a weakly coupled theory, such as electroweak theory. By 
comparison with Eq. pi) g"^ is analogous to T, and the formation of a QCD 
soliton condensate at T = is analogous to a high-T condensed-matter 
condensate such as a condensate of dislocations at crystal melting. 

2 Outline of the review 

In this review we focus on four QCD entropies at zero temperature, all of 
which are closely related. 

• Sohton entropy: Of greatest interest are the quantum solitons of an IR- 
effective QCD action (Sec. pi) that describes generation of a dynamical 
gluon mass. There are center vortices and nexuses in d = 2, 3, 4, having 
codimension 2, that are important for confinement and chiral symmetry 
breakdown (CSB). There are also the usual classical solitons, self-dual 
instantons and sphalerons, at short distances. 

• Condensate entropy: QCD has condensates such as (TrG'j^^,) and {qq). 
These are global descriptions of soliton condensates. 

• Entanglement entropy: Even a pure state, such as the QCD vacuum 
state, can lead to entropy if the Hilbert space is divided into two parts 
A and B with the pure state entangled between them, and a partial 
density matrix for part A is constructed by tracing over the observables 
inB. 

• World-line entropy and CSB: Nearly massless would-be Goldstone par- 
ticles, such as the pion, have Wilson loops that are highly branched and 
strongly fluctuating, because of entropy. Because the entropy gives neg- 
ative contributions to the action and mass of a particle it is critical for 
generating a massless pion when only (positive) confining effects are 
considered. Entropy also plays a critical role in generating the {qq) 
condensate, both of which accompany CSB. 

We have no space here to review the fundamentals of non-perturbative 
QCD that underpin these entropic studies; a few essential points are made 
in Sec. |3} The reader should consult recent books [U |2] for details of, e.g., 
the center vortex model and many references. 



Sec. |4] gives some general arguments for the entropy of condensates, not 
requiring detailed knowledge of their soliton content. Sec. |5] reviews the stan- 
dard lore for estimating the entropy of point-like solitons such as instantons, 
sphalerons, and d = 2 vortices. We go on to more complex estimations for 
center vortices in d = 3 that are analogous to closed polymers]^ A conden- 
sate of closed polymers maps onto a scalar field theory, and entropic effects 
contribute a negative term to the squared scalar mass. If there is entropy 
domination it becomes necessary to add a repulsive 0^ term for stability. 
In d = 4 center vortices are closed two-surfaces and described by an as-yet 
unknown string theory; we do not speculate on what such a theory might 
be like here. Sec. |6] reviews the canonical picture of confinement by center 
vortices, where the area law for a planar Wilson loop arises from an entropic 
calculation quite similar to that of the entropic rubber band. However, it 
is far from clear how to go on to more complex situations, such as a non- 
planar Wilson loop. In Sec. [7] we introduce a crucial ingredient missing in 
the techniques of Sec. |6] to show how for a non-planar Wilson loop the area 
law involves the area of a minimal spanning surface. The crucial ingredient 
is the difference between the piercing of any spanning surface by a center 
vortex and the actual topological linking of the vortex and the surface. 

The techniques of Sec. [7] have evident analogies to a genuine entropy in 
the QCD vacuum: Entanglement entropy. It arises even for a pure state 
as the residual entropy of this pure state when it is entangled between two 
subsystems A and B, and B (for example, the interior of a black hole) is 
unobserved. One of the most interesting properties of entanglement entropy 
is that for many systems, including black holes, entanglement entropy agrees 
with the holographic principle: It does not scale with the volume of the black 
hole, but with its horizon area (as does the Bekenstein-Hawking entropy). For 
QCD the pure state is the vacuum, and the unobserved regions will depend 
on the questions we are asking. In Sec. [8] we briefly review the fundamentals 
of entanglement entropy and go on to point out their relevance to the ideas 
of Sec. E 

Finally, in Sec. [9] we discuss the configurational entropy of massless quark 
world lines and their associated Wilson loops, which are important in how 
CSB and confinement are connected. For massless quarks this entropy con- 



^The same principles apply to nexuses in d = 4 (Nexuses are QCD magnetic monopoles 
whose flux tubes are embedded in center vortex surfaces. Nexuses disorient center vortex 
surfaces and give rise to topological charge when a nexus world line links to a vortex 
surface. 



tributes a significant negative term to the pion mass, making it possible to be 
a Goldstone boson with positive area-law and kinetic contributions. Entropic 
effects also drive the formation of a qq condensate. 

No discussion of QCD entropy would be complete without mentioning 
condensed matter analogs. There are lattice systems whose ground states 
are not unique, so their entropy does not vanish at T = 0. These include 
frustrated antiferromagnets,[3] water ice,|lj and spin ices. [SI E] They provide 
useful analogs to the QCD vacuum even though their idealized ground states 
are not unique. Disorder at T = arises because there is no energy cost 
to move from one realization of the system to another. For example, in 
a triangular antiferromagnetic lattice once any two spins in a unit cell are 
oppositely-aligned, the third spin may be anti-aligned with either, but not 
both, of the first two with the same energy for either anti-alignment. Spin 
ices actually have ferromagnetic interactions obeying the so-called ice rule, 
calling for two of the four spins in the tetrahedral unit cell to point in, 
while the other two point out. Even at zero temperature there is no unique 
way of arranging spins to obey the ice rule, hence disorder and entropy. It 
takes an extremely small energy to violate the ice rule {e.g., three in, one 
out) in one tetrahedron; a neighboring tetrahedron then flips a spin to obey 
the ice rule. This spin flip propagates along configurations that have disorder 
greatly resembling that of QCD, with quasi-particle magnetic monopoles and 
Dirac magnetic flux tubes joining them. Ordinarily a flux tube would cost a 
fair energy per unit length to form, but not for the spin ice. Although the 
physics is quite different, this is analogous to the balancing of entropy and 
action discussed in Sec. [3l 

3 The IR-effective action for QCD and cen- 
ter vortices 

In (i-dimensional non-perturbative QCD with gauge group SU{N), the dom- 
inant solitons appear in a gauge- invariant effective action [2] that summarizes 
the fundamental infrared quantum effect of dynamical gluon mass generation: [7] 

/ = ^ jd'x {\g-G + Tr'^iU-'D.Un, G-G^J2 ^U- (4) 

Throughout this paper we define the gauge potential and field as the canoni- 
cal ones multiplied by g. Here the trace is normalized to the sum over space- 



time and group indices, D^ is the covariant derivative, and path integrals are 
indicated by enclosing the differentials in parenthesesjj The traceless N x N 
matrix U transforms according to the adjoint of the gauge group under gauge 
transformations of the potential. The running coupling g"^ and mass m^ are 
each evaluated at zero momentum. 

This effective action has a number of solitons, among them center vortices 
with co-dimension 2. Vortex-like excitations, and the importance of entropy 
dominance, were invoked before center vortices [9] through simple Abelian 
lattice models. For non- Abelian gauge theories in d = 3,4, 't Hooft[lOj in- 
troduced center vortices that were soon thereafter interpreted [TT] as solitons 
of the effective action of Eq. dll). We will be explicit only for d = 3, where a 
simple center-vortex soliton has the form: 



^«(^j J> 



2'nQjtiabda f dzb[Ao{x - z) - A„(x - z)\. (5) 

Jcv 



The diagonal matrices Qj, j = 1 . . . A^ — 1 are such that exppyriQj] = 
exp[27ri/A^] is in the center of the group, and as before A^ is the free prop- 
agator of mass m. In other dimensions the integral is over a closed surface 
of codimension 2, and the epsilon symbol has d indices. This simplest of 
center vortices is essentially Abelian, but there are many ramifications, such 
as nexuses that disorient the vortex and contribute to fractional topological 
charge, that are essentially non- Abelian; [2] we have no space to study them 
here. Note that there is no singularity in the vortex on its defining contour 
CV ^ because short-distance singularities cancel between the two propaga- 
tors. In any dimension center vortices have finite transverse thickness 1/m 
coming from the massive propagator. The term in Aq, considered by itself, is 
a pure (singular) gauge transformation of long range; it alone is responsible 
for confinement. 

The dynamical gluon mass, absent at the Lagrangian level, is a conse- 
quence of condensates such as center vortices and nexuses. If there are no 
condensates there can be no gluon mass, since it has been shown [T^ that at 
large momentum this mass runs as: 

m\q-) ~ ^^ (6) 



^ We always work in Euclidean space at zero temperature, and suppress writing gauge- 
fixing and Faddeev-Popov terms. 



(just as the CSB running quark mass runs like {qq) /q^). And, of course, if 
there is no mass scale in the effective action there can be no solitons such as 
center vortices and nexuses. So the gluon mass, on which all non-perturbative 
effects depend, is entropy-driven. 

A major concern of this review is the area-law behavior of Wilson loops 
coming from their coupling to a condensate of randomly-distributed center 
vortices. We evaluate this for a single vortex and a Wilson loop F, all of whose 
lengths are large compared to 1/m. The massive propagator in Eq. ([s]) yields 
only perimeter terms and can be dropped, yielding: 

W = TrPe'frd-,M-) = exp[?^ hx^ (f dx,e,,fc5fcAo(x - z)] (7) 

^^ Jr Jcv 

where F is the Wilson loop contour and CV labels a given center vortex loop. 
The coefficient of 2'n'i/N in the integral on the right-hand side is a topological 
quantity, the Gauss linking number]^ so for a single vortex 

W = e^-^'^^l^. (8) 

The integral over the loop contour can be rewritten with Stokes' theorem to 
become an intersection number of the center vortex with any surface that 
spans the Wilson loop. This raises interesting questions about which area 
appears in the expectation value (VT), which will be taken up in Sec. [7J 

The next question is how to describe the entropy of a condensate of center 
vortices or other solitons. 



4 General arguments for condensate entropy 

Is it possible to say anything about the entropy of a condensate, without even 
knowing what the condensate is made of? Perhaps surprisingly, the answer 
is yes. Consider first an often-discussed candidate for a "condensate" that 
has zero entropy: A covariantly-constant chromomagnetic field in d = 4. Use 
the notation 

Q{x) = \g-G (9) 

where h = lliV/487r^ is the no-quark one-loop coefficient in the beta-function 
(3{g) = —bg^ for SU{N). The one-loop action can be found [T3] by the stan- 



■^This is true in all dimensions. 



dard techniques of Schwinger, and yields for the real part: 



r{e) = jdxein{—}. (lo) 

This is an interesting expression because it has a minimum of negative value 
— J{9) at 9 = (9), suggesting the dominance of entropy over (positive) ac- 
tion. Unfortunately, it has been pointed out[Tl] that, just as in the Schwinger 
calculation of vacuum decay for a constant electric field in QED, one should 
add to this action an imaginary part, signaling instability of a constant chro- 
momagnetic field in QCD. Nielsen and 01esen|Ll4j showed that the instability 
could be removed with sufficient space-time inhomogeneity, and correctly in- 
terpreted this as a signal for formation of high-entropy disordered domains. 
But one need not assume a covariantly-constant field to arrive at F of 



Eq. (10). We can find this action, with no imaginary part, using methods 
based on the one-loop renormalization group. p^ [T6] Of particular interest 
is an infinite set of relationsflG] among zero-momentum components of 9: 



N 



J ... jifl dx,] {T9{x,) . . . 9{xr,)9{0)),onn = {9(0)) (11) 

found by repeatedly applying g^d/dg to the vacuum energy density e^ac = 
(1/4) (T^) as defined by the partition function: 

Z(0) = e^-''^^^"- = f{dA^a)e'^^^'^ (12) 

from which it follows that 

9'^e,ac = ^{G-G). (13) 

Because eyac is renormalization-group invariant of dimension 4 it satisfies 

if^g- + m-g-^vac = (4 + mg-)^vac = 0, (14) 



which, with Eq. (13), yields the trace anomaly. Repeated differentiation of 



Eq. (12) then yields the sum rules of Eq. (11) 



At this point we have made no use of the putative one-loop effective action 



of Eq. (10). But standard Legendre transform methods show [17] that r(^) 



in Eq. (10), with no imaginary part, is precisely the effective action which 



gives rise to these sum rules. Moreover, the generating functional found by 



exponentiating Eq. (10) is[T8] a Poisson distribution for 9, if this quantity is 
thought of as the sum of a large number of independent random terms with 
average value (9), so that 9 = J^^i- This suggests strongly that the {G ■ G) 
condensate has a large entropy, arising from lumps of condensate with a finite 
correlation length. 

These considerations depend critically on using the beta-function at one 
loop only. But in (i = 3 (where there is no beta-function) the infinite sum 



rules analogous to those of Eq. (11) and the corresponding effective action 
are exact, given the reasonable hypothesis that d = ?> QCD has only one 
mass scale, g^ itself. [18] In c? = 3 we define 9 as: 

9{x) = ^^G-G, (15) 

an operator with mass dimension 3. Because g^ has mass dimension 1, e^ac 
necessarily scales like g^ in the absence of any other mass scale. This might 
seem impossible since in perturbation theory high-order graphs have high 
powers of g. But all graphs of eighth or higher order are UV finite; their IR 
divergences are removed by a gluon mass m ~ f?^, so ordinary dimensional 
counting works. Eqns. (T2p!3) continue to hold in d = 3, and combined with 



Kg for some pure number K yields: 

1 



vac- ^{0). (16) 



Since (9) is positive, the vacuum energy is negative and hence entropy- 
dominated]^ As in (i = 4, continued application of g^d/dg to e^ac yields 
an infinite set of sum rules: 

j ... jifl dx,] {T9{X,) . . . 9{xM)9{0)),onn = ^^^^^' (^(0)). (17) 

Again, standard Legendre transform techniques yield the effective action: 

T{9)= U\{9-U'"\9Y/^). (18) 

"'This equation is actually a simple consequence [18^ of Lorentz invariance and isotropy 
of the vacuum in d = 3. 



This has its mininiuni value of — (l/3)(0), as required. 

There is[T9] a complex scalar field theory with a {ipl'^ interaction in which 



loop corrections yield a vacuum energy of precisely the form of Eq. (18), 
with the identification 6 ~ \ip\^- The peculiar 3/4-power term comes from 
loop corrections yielding (specific to d = 3) \ip\^ terms that, it turns out, 
are equivalent to a negative squared mass. Just such a theory is needed to 



describe a d = 3 center- vortex condensate, as we review in Sec. 5.2 



5 Entropy of solitons 

5.1 Entropy of instantons 

The simplest expression of entropy is for point-like solitons, such as instan- 
tons. Strictly classical instantons, the solitons of the usual (m = 0) action, 
can have any size. But in the IR-effective action these are approximately 
solitons only if their size p is very small in the sense mp <C 1. Otherwise, 
the solitons analogous to instantons (which are no longer self-dual) have size 
limited by 1/m. Similarly, the barrier height represented by sphalerons can- 
not vanish, but is of order m at least. For simplicity, and because we are 
interested only in IR effects, we will proceed as if instantons were of fixed 
size and ignore the size collective coordinate. 

Set the currents J to zero in the generating functional equation ([3]), since 
they do not contribute to the entropy, and add a chemical potential term 
exp[pN] in the summand over N. This generating functional at zero current 
is just the grand partition function at zero temperature. The only integral 
over collective coordinates that we need indicate explicitly yields a factor 
V4^/Vc for every instanton, where Vc is a. finite volume (depending on g'^)to 
be calculated from the determinants. Then, using Stirling's formula and 
replacing the sum over N with evaluation of the summand at A^ = iV ~ V4 
we find 

\nZ = N[ln{^J-I^ + p] (19) 

This is equivalent to the standard thermodynamic result InZ = P[—U + 
TS + pN] for free particles with internal energy, and so we can identify the 
entropy per instanton as the logarithmic term in the brackets, and the free 
internal energy per particle, in units of temperature T, with the classical 
action L. 
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It was once argued that instantons (and by extension, other sohtons) 
could not contribute much at large Nc for color gauge group SU{Nc), because 
of 't Hooft scaling g'^ ~ 1/Nc and the action factor is suppressed exponentially 
in Nc- But in fact the leading term in the collective coordinate integral 
exactly cancels this action suppression. [20] What happens at non- leading 
orders is unknown, but there is no known reason to believe that instantons are 
suppressed at large Nc- The same leading-order exact cancellation between 
action and entropy also takes place for center vortices, at least in d = 2, 
where they are point particles. [8] 

More elaborate versions of the above argument were used early on to 
argue for confinement by instanton condensates, but these ultimately failed 
for a number of reasons that we will not go into here. Perhaps it should be 
noted that recent lattice simulations suggest that dilute instantons and anti- 
instantons do not constitute a large fraction of the QCD vacuum configurations. [21 

5.2 Entropy of center vortices 

Dealing with the configurational entropy of center vortices is conceptually 
straightforward in d = 3, where a condensate of oriented closed strings can 
be mapped onto a complex scalar field theory (a trick long used by polymer 
physicists). The idea[22] is that the partition function for a condensate of 
these strings is exactly the same as the partition function for the free complex 
scalar field theory, once one recognizes that the propagators in the vacuum 
graphs of the scalar field theory can be written in Schwinger-Feynman proper- 
time form as an integral over closed strings. First, the standard partition 
function for a free complex scalar theory of mass M is 

Z = / (dV'dV')e-^« =exp[TrlnAM] (20) 

Jo = /d^x |(9^|2 + MVr 



where Am{x, y) is the free propagator, and the Tr is over all Euclidean space- 
time (and possibly other discrete degrees of freedom). Up to a normalization 
factor the trace of the logarithm of a scalar propagator (in any dimension) 
has the proper-time representation as a path integral over closed world lines: 

Trln[AM] = y ^j{dz)exp{-j dT[-z' + —]}. (21) 

11 



Now note that in d = 3 the probabihty P{x, x'; t) that a random string made 
up of t steps, of average length a, will go from x to x' has the continuum 
random-walk path-integral form 

Pix,x';t) = j (D^)exp[-y dr— ] (22) 

(up to a normalization) where dr ^ a. Here a represents the correlation 
length for the string to change directions by, say, tt/2. 

The crucial step in estimating the entropy is to find the number of ways 
to go from X to x' with a string of length t. Imagine that the string lies on 
the links of a cubic lattic of side a. The string grows by addition of a new 
link at one end, which can be done in 2(i — 1 = 5 ways when backtracking 
is disallowed. Then the number of ways to go from x to x' is approximately 
P{x,x';t) times exp[51nt]. We weight each link of the string by a factor 
exp[— /i], or exp[—fit\ for the whole string; for the field theory propagator 



this represents the mass factor in Eq. (21) 



The weighted total number of distinguishable ways a single oriented closed 
string going from a; to x can be formed is roughly 

5*/ dT-P{x,x;T)e-''' = Zi (23) 

Jo 'T 

which defines the one-string partition function Zi. The factor 1/r corrects 
for double counting: Divide by the length of the string, because any point 
on a given string serves to label a single configuration. The grand partition 
function is, as usual, 

Zstring = ^ ^^ = exp Zi. (24) 

We are almost done. Change variables via ds = (a/3)dr, and then define 

M^ = [4][/i-ln5]. (25) 

The resulting grand partition function Zstring is precisely the free scalar par- 



tition function of Eq. (20). Observe that entropy can dominate over the 
standard mass term, in which case M^ is negative, and the path integrals 
defining the scalar propagators make no sense. This is, of course, just the way 
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the mass term behaves in spontaneously-broken theories. Such theories have 
to be stabihzed by a \tp\'^ interaction, and then there is a phase of the theory 
in which (|^|) is a non- vanishing constant throughout space-time, signahng 
a condensate. In our case we interpret the l?/'!^ interaction as mimicking the 
mutual avoidance of different strings, since the overlap action of two string 
trying to occupy the same space-time amounts to a repulsion. The conden- 
sate forms because the entropy per unit vortex length dominates the action 
per unit length. Then vortices tend to become very long; in fact,^] a finite 
fraction of them have infinite length. 

From the explicit form of the d = 3 center vortex in Eq. ([s]) the expecta- 
tion value (W) of a Wilson loop T is equivalent to coupling a gauge potential 
to the field ip. For the simple case of SU{2), where there is only one ma- 
trix Qj with entries ±1/2, the vortex is effectively Abelian with Qj = 1/2 a 
number and not a matrix, because every configuration is accompanied by an 
equal contribution from its complex conjugate. This vortex is: 

Vi{y) = 11 j) dxjeijkdk[^Q{x - z) - Am(2; - z)]. (26) 

If one is only interested in the area-law part of the Wilson loop, the A^ 
propagator of Eq. ([s]) can be dropped. 

This model's implications for confinement and for Chern-Simons fluctu- 
ations have been explored|19j. beginning with a complex scalar field theory 
with \iIj\'^ coupling. Quantum corrections yield the results of Sec. E^for en- 
tropy generation in rf = 3. These considerations have been generalized [23] 
to non- Abelian gauge theories beginning from the center vortex solitons of 
Eq. ([5]) and adding nexuses. The partition function of the IR limit of the 
gauge theory has the form of the string partition function of Eqs. ( 23p4 ) 
after a simple redefinition of variables. 

Four dimensions is much more complicated, since the center-vortex par- 
tition function naively maps onto a string theory in a non-critical dimension. 
It will be essential to modify the action of this string theory with curvature- 
dependent terms accounting for stiffness and bending resistance of the vortex 
surfaces, coming from increases in the fundamental gauge-theory action when 
a vortex surface is bent. However, these complications are not particularly 
important for understanding the basics of confinement, which we present in 
the d = ?> context. 
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Figure 1: A rectangular Wilson loop whose flat spanning surface is divided 
into A^ squares whose length is the correlation length A. Either zero or one 
vortex can pierce any square (black dots indicate a piercing 



6 Entropy and confinement 

The Wilson-loop area law is actually entropic. As remarked in Sec. [3} a single 
center vortex contributes a phase factor 



W = Pexp[i (bdxiAi 



exp[ 



,2'KiLi 



N 



(27) 



where Lk is the Gauss linking number of the vortex tube and the curve F. 
The phase factor for many vortices is the product of the phase factors of 
the individual vortices, which do not interact. For simplicity consider only 
SU{2), where the fundamental-representation Wilson loop for a single vortex 
is (— l)^*. Even Gauss link numbers contribute nothing to confinement; only 
odd links count. 

To find[Tl 12] the area law of (W) choose a simple loop, consisting of a 
rectangle spanned by a flat surface (see Fig. II]). The spanning surface is 
divided into A^ squares whose side is the correlation length A, so that the area 
A of the square is ^ = A^A^. This length is (roughly) defined by saying that 
no more than one vortex can pierce, or intersect, such a square, so a square 
can either be singly-occupied (with probability p) or empty (with probability 
1 — p). Because we use the DGA, we take it that p ^ 1. 

The black dots indicate where a vortex has pierced the Wilson-loop sur- 
face; we assume such piercings are random and independent of one another. 
Although a piercing is not necessarily the same as a linkage (that is, L^ is 
odd) to the loop, because it simplifies things so much we assume that ev- 
ery piercing indicates a link number of unity. In that case, the expectation 
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value {W) is just finding the probability P/v( J) that there are J linkages, so 
defined, and calculating the expectation value of (—1)'^ with this probability. 
Here is the first place where entropy enters, and it enters in precisely the 
same way as it does for the entropic rubber band. By simple combinatorics 
the probability Pn{J) is 

the same as for an anisotropic rubber band, where the probability p of a left- 
pointing link is less than that of a right-pointing link. The factorial factors 
are just the number of ways of distributing J things in N slots. To relate this 
probability density to an entropy density, note that there are a total of 2^ 
states, with an entropy of A^ In 2. The number of states for a given occupancy 
number J is the product 2^Pn{J), whose logarithm defines a partial entropy 
at fixed J. In expectation values this 2^ factor cancels with the same factor 
in the partition function. 
We easily calculate 

(W) = E i-^yPN{J) = (1 - 2p)^ = e-^--^ (29) 

J 

where 



Kf = r^ ^ — = 2p2 (30) 



ln(l - 2p) _ 2p 

and p2 is the areal density of pierce points. [Because the codimension of a 
vortex is 2, vortices are characterized in any dimension by an areal density of 
pierce points, and a confining Wilson loop is an area law in any dimension.] 
Note that the areal density is necessarily the same anywhere in spacetime. 

However enticing the simplicity of this "derivation" of an area law may 
be, it has a serious fiaw, which when resolved will reduce the entropy: The 
piercings are correlated. It is much more complicated to deal with this issue, 
which as it turns out does not interfere with the basic fact that confinement 
means an area law. [23] Furthermore, the derivation does not explain the 
"obvious" point that the surface whose area appears in the area law is a 
minimal spanning surface. How does one show that for a non-planar contour? 
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Figure 2: Two surfaces spanning a non-planar Wilson loop (heavy line): A 
minimal surface M, and a much larger surface L. How to tell which area, M 
or L or either, to use in the area law? 



7 The difference between piercing and link- 
ing 

One should not conflate piercing with linking. There are critical circum- 
stances in which a vortex necessarily pierces a spanning surface an even 
number of times, in which case it makes no contribution to the area law. 
Nonetheless even-numbered piercings are important, because they render in- 
active portions of the spanning surface that might otherwise be linked. This 
inactivation changes the entropy crucially for non-minimal spanning surfaces. 
(There is also a certain amount of inactivation even for a minimal surface, 
but it is a quantitatively non-leading effect. p^) 

Draw a generic (i.e., non-planar) large Wilson loop such that it has a 
unique minimal spanning surface M, and then erect a different spanning 
surface L with necessarily greater area so as to form a closed surface on 
which the Wilson loop lies (see Fig. d2])). The areas of these surfaces are 
•Amj-^l- (For reference to Sec. [SJon entanglement entropy, we now have a 
closed surface M-l-L dividing three-space into two parts A and B, with B the 
region interior to the surface M+L.) Which area are we to use in the Wilson 
loop? Since the contribution of a given vortex to (W) is a purely topological 
signed intersection number, it is the same, and more important it is also the 
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same (mod 2), for either spanning surface M or L. 

Is there a unique answer for the area law coming from vortices? The 
answer is yes, and explains why the area law involves a minimal spanning 
surface. The demonstration only uses knowledge of what happens on the 
surfaces M and L. Only vortices entering the volume B by piercing M and 
exiting by piercing L (or vice versa) give a linkage with the Wilson loop. 
There are, by assumption, N^ = P2AM vortices piercing M. By the same 
assumption Nl = P2-Al vortices pierce L, and Nl > Nm- In fact, Nm is 
the least number of vortices (up to small fluctuations) that can pierce any 
spanning surface. It should therefore be clear that there are {N^ — Nm)I'2 
vortices piercing surface L that do not pierce surface M at all; they exit 
through surface L again. Their only role is to render inactive enough vortex 
sites so that exactly Nm are really linked. So in using surface L to calculate 



the area law by invoking Eq. (29 ), the number of vortices entering this formula 
is Nmi not Nl, so we find the unique minimal area law. A study |23] of the 
question glossed over here, of the quantitative difference between the number 
of vortices piercing a flat surface and the number linked to it, shows that for 
dilute vortices the string tension is slightly reduced from its canonical value 
of 2p2. 

Here is an explicit computation of the area law for a simple non-minimal 
surface, showing the importance of distinguishing piercing and linking. Change 
the spanning surface M to be slightly non-planar by excising one interior 
square, and erecting a cube of side one unit over this square; the bottom 
square of the cube is missing. We call this cube the bump. To generalize the 



entropic calculation of Eq. (28) use the identity: 

N 

j](Ai+5i)= j]5,+^A,j]5,+j]AA, n^-^+'-'+n^- (31) 

Suppose that every square in the lattice of Fig. [T] has a different probability 
Pi of piercing and a probability 1 — pj of no piercing. If we identify Ai = 



Pi,Bi = 1 — Pi, the right-hand side of Eq. (31) decomposes the total (unit) 
probability into terms with no, one, two,. . . piercings. To find (W) for this 
lattice, just set Ai = —pi, Bi = 1 — pi. The result is: 

N 

{W)=l[{l~2p,)^e-^^P^ (32) 



which, of course, agrees with Eq. ( |29[ ) when pi — > p. For the lattice with the 
bump, there are A^ — 1 squares with probability pi = p, and 5 squares with 
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probability p/5. This follows because if a vortex (which is now allowed to 
form a three-dimensional closed path along the dual lattice) pierces any of 
the 5 squares in the bump, there is only one way to have an odd link number, 
and that is to go out the empty bottom square. All other exits have even 



(in fact, zero) linking number. Using these results in Eq. (32) yields exactly 



the result in Eq. (29) for the fiat lattice without any bump. If we had just 
counted piercings, without asking whether they amounted to linkings, this 
would certainly not be so. It is more tedious to do these calculations for 
various other bumps (two adjoining cubes,...), but one always gets the same 
result: 1) Any two surfaces spanning the flat Wilson loop give the same area 
law; 2) This area law is that of the minimal surface. 

Our use above of the closed surface M+L that bounds a three-volume 
where we know nothing about what happens to the center vortices in the 
interior strongly suggests an analogy with entanglement entropy, [251 IMl [271 
|28| |29| 1301 l3H 132] as pointed out in Sec. [s] immediately below. As before, 
the calculation of an actual area law is the same as calculating the entropy 
of vortex link points on these surfaces. 

In view of the analytic difficulties in distinguishing piercing from linking, 
it would be worthwhile to do lattice simulations by constructing ensembles of 
random mutual- and self-avoiding vortices on a lattice dual to the lattice on 
which Wilson loops and surfaces five, and calculating expectation values of 



the type in Eq. (29) numerically for a variety of non-planar Wilson loops, 
to establish that the area law holds with the minimal area and a universal 
(loop-independent) string tension. One can then also study fluctuations in 
the area of the area law, to establish such effects as the Liischer term. These 
simulations will be far easier than simulations of full QCD. 

8 Entanglement entropy 

Usually one thinks of the entropy of a pure state \ip) of a system as being 
zero. However, if the system can be decomposed into two subsystems A, B 
whose sum is the whole system, and one of the systems is unobservable for 
some reason (for example, B is the region inside the horizon of a black hole), 
the partial entropy of the remaining subsystem will not vanish even for that 
pure state. This can happen when the pure system state is an entangled 
state of the two subsystems. We define a density matrix and entropy for the 



A subsystem in the usual way: 

Pa = TiB Pa+b = TiB I ^)(^ I, S'a = - Tr a Pa In pa- (33) 

Suppose A and B are two-level systems, and we denote their eigenfunctions 
in the usual way as |+) and |— ). We calculate the entanglement entropy of 
the pure but entangled state 

\ip) = cosa\ H — ) + sinal — h) (34) 

(state of subsystem A written first in the ket) from 

Pa = cos^ a|+)(+| + sin^ a|— )(— I, (35) 

Sa = — cos^ a In cos^ a — sin^ a In sin^ a. 

This vanishes if cos a or sin a vanishes, when the state is not entangled. 

Entanglement entropy has many applications in condensed matter theory. [271 
I2H1 E2] field theory, [23 [il ED W\ and string theory. [23 A standard example 
is a nearly-static black hole in d = 3 -|- 1 space-time whose horizon divides 
three-space into two regions A (observable) and B (unobservable). As is well- 
known, the entropy of the black hole does not scale with the volume inside 
the horizon, but with its area. 

We have already seen such a feature arise in QCD in Sec. [7[ where the 
closed surface M + L of Fig. |2] divides three-space into outside (A) and inside 
(B) regions. The entanglement entropy scales with the area of region B. The 
observable (that is, in region A) center vortices of the vacuum are closed 
loops, but those closed loops entering and leaving the surface of region B 
appear to be open strings in A, terminating on this surface. Such loops 
are, of course, entangled with their continuations in region B. Just from 
observations in region A, one sees the open strings and their endpoints on 
the surface, very much the same as if the surface were a brane with the 
chromomagnetic flux tubes of center vortices terminating on itjj We have 
already calculated the entropy associated with the random positions of these 
endpoints, which behave for all practical purposes like a gas of point particles 
on the surface, with the usual vortex areal density p2, and found a unique 
minimal-surface area law. 



"■The identical argument has been given [31j but with the postulation of a disordered 
condensate of chromoelectric rather than the chromomagnetic flux tubes of the center- 
vortex condensate. 
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Another example [5^ of entanglement entropy (not termed as such) is 
based on a d = 3 analog of the original Liischer bag[35] that divides four- 
space into an outside A and an inside B. In this example the Chern-Simons 
term plays the role in the three-dimensional bag that the G^yG^y density does 
in the Liischer bag. We interpret this bag as region B, with an unobserved 
Chern-Simons density. The volume integral of this density can be expressed 
as a surface integral that expresses the holographic entanglement entropy. 
Entanglement entropy from long-range topological disorder is also important 
in condensed matter physics. [28j 



9 Entropy of Wilson loops: Chiral symmetry 
breaking 

The idea of chiral symmetry breaking (CSB) goes back to the fifties, well 
before QCD, yet it is still an active field of research. Lattice simulations have 
helped, but not singled out a definitive mechanism. For SU{2) they show [36 j 
that center vortices not only confine quarks, but are also critical for CSB: 
When they are artificially removed, both confinement and CSB are gone. 
This is consistent with earlier continuum arguments j37] that confinement 
is necessary and sufficient for CSB. There is also lattice evidence that the 
deconfinement transition temperature is rather close to the CSB restoration 
temperature for quarks. 

For decades people have embraced the idea that one-gluon exchange is 
strong enough to produce CSB without confinement; we examine the opposite 
case, where confinement is needed. Note, however, that a recent study [3S] 
concludes that one-gluon exchange is enough for CSB, based on a study of 
the pinch technique [2] Schwinger-Dyson equations based on lattice results 
and including a gluon mass, which tends to weaken one-gluon effects. 

9.1 The importance of entropy for the pion 

Up to now the Wilson loops that entered our study of confinement were very 
large and prescribed beforehand. Since they were fixed, they had no entropy 
at all. The physical interpretation of such large smooth Wilson loops is that 

1 /9 

they refer to very massive particles, with mass Ai ^ Kp [Kp is the string 
tension). Such particles can be far apart (of order M./ Kp) and hardly feel 
the confining force. 
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We now want to study dynamical Wilson loops in the opposite limit, of 
massless quarks that make up the massless Goldstone pion, a q'j^q bound 
state. [371139] We work in the quenched limit and omit quark determinants in 
path integrals. In the quenched limit the propagator of the pion field 7i{x) 
is: 

(T7r(a;)7r(0)) ~ /"(dA)e-^(^)Tr[5^(x, 0)75^a(0,x)75] (36) 



where Sa is the quark propagator in the gauge potential A; it is an integral 



over paths, in the same spirit as Eq. (21). The trace creates a Wilson loop 



starting at 0, passing through x, and returning to 0. Between these endpoints, 
the Wilson loop fluctuates, and if the mass parameter in the propagator is 

1/2 

small enough compared to Kp these fluctuations are large and have large 
entropy that will make important negative contributions to the mass of the 
pion. Such negative contributions are essential, if one-gluon exchange is (as 
we assume) too weak to produce CSB by itself. 

The distance r between q and the q in the pion cannot be large; if it were, 
the confinement energy Kpr would be considerably larger than such scales 
as the running quark mass at zero momentum M(0). On the other hand, the 
qjbq Wilson loop can propagate a long distance, because the pion it represents 
is massless. So schematically we can think of the pion Wilson loop as very 



long and thin, and apply the ideas used earlier in Sec. |5.2| for the entropy 
of center vortices. Based on a simple picture [3Z] of the confining dynamics 
of a massless qq bound state, a good guess for the inverse correlation length 
between bends of the pion world line, which is essentially the entropy per 
unit length, is the inverse of the distance at which the confining energy is 
comparable to the quark mass, or Kf/M{0). This is quite comparable to the 
action per unit length, and should make an important negative contribution 
to the pion mass. It is, in fact, indispensable if we ignore one-gluon effects, 
because the other contributions — the Kpr confinement energy and the kinetic 
energy — are positive. 

Another important entropic effect [STf [3U] is the sprouting of long thin 
branches off the main pion trunk, representing the qq condensate that sig- 
nals CSB. Such a condensate can only form [37] if the g or g world lines have 
segments going backwards in time in Minkowski space, which is highly im- 
probable for heavy quarks but not for massless ones. The final picture of the 
pionic Wilson loop is that of a highly ramified and convoluted one with no 
large spatial separation r between the q and q. Such a Wilson loop is drasti- 
cally different from the ones invoked earlier in Figs. QJ pi and seems to be|40j 
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in the equivalence class of branched polymers, with Haussdorf dimension 
4.|ll] 

It is not easy to deal directly with dynamical Wilson loops, highly- 
ramified or not (see an example ^5] of such a CSB study, with further ref- 
erences). Far better would be some form of Schwinger-Dyson equation that 
captures entropic effects more or less automatically. Next in Sec. 9.2 we, as 
many others do, use a simplified approximation to an area law expressed as 
an effective vectorlike propagator ~ 1/k'^ coupled to the quarks. When the 
quark worldline is closed, as it is in a Wilson loop, there is an Abelian gauge 
invariance that has the remarkable effect of automatically regularizing the 
IR divergence of this propagator with a physical cutoff. Moreover, it turns 
out that this physical regularization gives rises to a negative contribution to 
the pion energy that allows for CSB. 

9.2 An area law and entropy 

In principle, the minimal-surface area law is expressible as a functional of 
the Wilson loop contour. Only in d = 2 is this simple, because the minimal 
surface is flat: 

Ar = (p dzi (p dzj -r-^ij In \z — z\. (37) 

But in c? = 3, 4 the contour form|12] of the minimal area is far from elemen- 
tary, and difficult if not impossible to implement (see, however, studies |15| Hi] 
of this problem). So we continue to use the simple form above, adjust- 
ing the overall factor as needed in c? = 3, 4. In dimension d the Fourier 
transform of the logarithmic propagator behaves like k~'^, and leads to IR 
singularities when used naively for such applications as CSB. For example, 
the conventionally-defined static potential from this propagator can only be 
defined by supplying an infrared cutoff, and then this potential rises like r as 



needed. A closely-related issue is that Eq. (37) has an Abelian gauge invari- 
ance that has nothing to do with color gauge invariance. If terms ~ didj are 
added to the delta-function they give no contribution, because the contour 
is closed. Among other things, this means that we need not specify the scale 
of 1 2; — 2;' I in the logarithm because changing it is a gauge transformation. 



Without this gauge invariance there is no definite meaning to Eq. (37). Note 
that this potential, mimicking an area law, is always to be used to first order 
only in graphical applications. 
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The standard gap equation for CSB refers to an open quark line: 
where the dots stand for omitted gauge- dependent longitudinal terms. Pos- 



tulating the d = A version of Eq. (37) for the effective propagator 



fl.„(*)..-VB.«(*); fle„(*) = ^ (39) 

in this gap equation violates the Abelian gauge invariance and requires an 
IR cutoff if CSB is to occur, because if M(0) ^ there is an IR divergence 
in the integral, unless a special Abelian gauge is chosen. There is also, of 
course, the problem that this gap equation refers to an amplitude that is 
not color-gauge-invariant either. Whatever the Abelian or color gauge, one 
cannot be certain that the right physics is captured in those gauges. 

QCD physics should be restricted to studies of closed world lines, corre- 
sponding to a physically-realizable qq pair such as describe the pion or the 
CSB condensate, and possessing both Abelian and color gauge invariance. 
Nevertheless it should be possible f39] to extract an effective gap equation 
from a closed-loopjj amplitude involving not only the massless quark under- 
going CSB but also a hypothetical very heavy spectator quark. Imposing 
these two gauge invariances has a perhaps unexpected bonus: It automati- 
cally provides for a resolution of the IR divergence. 

First we take up the Abelian gauge invariance. The sum of all closed- 
loop massive-particle amplitudes to a given order (in our case, first order) in 



the propagator of Eq. (39) has this invariance. Call the sum of closed- loop 



graphs F, and construct the functional derivative 

By Abelian gauge invariance this is conserved: kaGapik) = 0; a direct calcu- 
lation verifies this. A standard kinematic argument shows that if (as CSB de- 
mands) there are no massless lines in F, Gapik) must vanish at least quadrat- 
ically at small k. When F is recovered by multiplying Gaii{k) by the effective 

^The gap equation is approximately what one would get by using the gauge-invariant 
pinch technique quark propagator j2j. 
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propagator and integrating over k, all IR divergences are removed, and re- 
placed by the lowest mass scale M = M(0) of the lines in F, where M(0) is 
the CSB-generated running quark mass at zero momentum. The would-be 
divergences are in effect "regulated" by this mass. A similar screening of an 
area law was found [46j in Coulomb gauge. This sort of screening should not 
be confused with many applications [17| of a screened area law having noth- 
ing to do with closed loops or entropy, but simply used as a device to mimic 
string breaking. An effective wavj39] of representing this physical regulation 
is to replace the original confining propagator by: 

D^ffik),^^5,^D^ffiky, ^e//(/^) = ^^^^^; (41) 

we expect m to scale with M(0). The static potential V{r) for this propagator 
is now finite: 

V{r) = -—e-'^' = -— + KFr + ... (42) 

m m 

where the omitted terms should be irrelevant because large separations r are 

never reached. 

The interesting thing about this equation is that the positive term 2M has 

been turned into a negative term! We suggest [39] that it represents entropy, 

similar to way in which the entropy of polymers or center vortices is effectively 



a negative contribution to the mass, as shown in Eq. (25). The regulator term 
is negative, and we take this as a sign of entropic forces at work. To find 
a useful estimate of m, consider the Hamiltonian H = 2M + V{r) + . . . 
for a color-singlet closed loop representing a qq bound state. Long ago it 



was shown [37] that with the original confining propagator of Eq. (39) both 
the mass term 2M and the static potential V{r) had IR divergences, their 
sum does not. But only recently [39] was the finite term remaining after 
cancellation of the divergences evaluted: 

2M + V{r) = Kpr -^. (43) 

Once again, an unexpected negative term appears. Equate the negative term 



in Eq. (43) to that of the regulated static potential of Eq. (42) and find: 



vrM 
m = ^ M. (44) 

Now we have tamed the IR divergences of the effective confining propagator 
with the physical quantity M. 
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These negative terms are quite essential for CSB with only confining 
forces, since the pion must have zero mass. A schematic Hamiltonian for the 
pion is: 

H^^2M + KFr + -^- + KFr-^^. (45) 

r r TiM 

Here we estimate the kinetic energies by 2/r. Without the negative (entropic) 
term a massless pion would be impossible, but with it there can be a massless 
bound state. A simple minimization of H.,; on r yields a zero eigenvalue at 

M2 = QKf/S. 

Now we can construct |39] a gap equation that is free of IR divergences, 
simply by using for Df,ff{k)^y the regulated propagator of Eq. (41), with 



m ^ M. The solutions automatically break CSB, and yield M^ k. O.SKp/ir, 
OT M ^ 250 MeV (after the inclusion of about at 10% increase from one- 
gluon exchange). Note that attempting to go smoothly by variation of some 
parameter from CSB to chiral symmetry restoration is impossible; the only 
dimensionless parameter is the effective coupling Kp/M'^, which shows un- 
bounded growth as one tries to take M to zero. This is quite unlike CSB 
with one-gluon exchange, where CSB is turned off by reducing the gluonic 
coupling. 

These results for QCD have been verified [i8] and extended|l9] to a [qq^ 
model for CSB in technicolor models. 
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